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Abstract. Infrared (IR) dust emission from galaxies is 
frequently used as an indicator of star formation rate 
(SFR). However, the effect of the dust-to-gas ratio (i.e., 
amount of the dust) on the conversion law from IR lu- 
minosity to SFR has not so far been considered. Then, 
in this paper, we present a convenient analytical formula 
including this effect. In order to obtain the dependence 
on the dust-to-gas ratio, we extend the formula derived 
in our previous paper, in which a theoretical formula con- 
verting IR luminosity to SFR was derived. That formula 
was expressed as SFR/(Mq yr^^) = {3.3 x 10-^°{1 - 
7y)/(0.4 - 0.2/ -t- 0.6e)}(LiR/LQ), where / is the fraction 
of ionizing photons absorbed by hydrogen, e is the effi- 
ciency of dust absorption for nonionizing photons, rj is the 
cirrus fraction of observed dust luminosity, and Ljr is the 
observed luminosity of dust emission in the 8-1000-/im 
range. Our formula explains the IR excess of the Galaxy 
and the Large Magellanic Cloud. In the current paper, es- 
pecially, we present the metallicity dependence of our con- 
version law between SFR and Lir. This is possible since 
both / and e can be estimated via the dust-to-gas ratio, 
which is related to metallicity. We have confirmed that 
the relation between the metallicity and the dust-to-gas 
ratio is applied to both giant and dwarf galaxies. Finally, 
we apply the result to the cosmic star formation history. 
We find that the comoving SFR at z ^ 3 calculated from 
previous empirical formulae is underestimated by a factor 
of 4-5. 
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1. Introduction 

When we want to know the entire evolutionary history of 
a galaxy, its color and metallicity are important diagnostic 
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quantities (e.g., Tinsley 19721 ). The evolution of the color 
and the metallicity results from the superposition of suc- 
cessive star formation during the lifetime of the galaxy. 
Thus, if we wish to reveal the evolution of galaxies, we 
always need to estimate the star formation rate (SFR) on 
a galaxy-wide scale. There are many methods to estimate 
the SFR of galaxies from observational quantities (Ken- 
nicutt 1998a). In this paper, we are especially interested 



in the conversion formula from infrared (IR) luminosity of 
the galaxies to their present SFR. The IR luminosity orig- 
inates from dust that is heated by stellar radiation field. 



Kennicutt (19981:) estimated SF R fro m IR luminosity for 
starburst galaxies. Buat & Xu ( |l996|) adopted an empir- 
ical approach by utilizing the observed relation between 
the ultraviolet (UV) luminosity and the IR luminosity. 



Recently, Inoue et al. (2000a, hereafter IHKGO) have 
derived a theoretical conversion formula from IR lumi- 
nosity to the SFR by developing a standard model of H 



II regions by Petrosian et al. (1972). The present paper 
examines how metallicity (or dust-to-gas ratio) affects the 
conversion factor of IHKOO. It is possible because the effect 
of the metallicity is well parameterized in their formula. 
Here, we present it as 



SFR 



Mq yr- 



3.3 X 10-i"(l-?7) £iR 
G.4-0.2/ + G.6e Lq 



(1) 



where / is the fraction of ionizing photons absorbed by 
hydrogen, e is the efficiency of dust absorption for non- 
ionizing photons, LiR is the observed luminosity of dust 
in the wavelength range of 8-1000 /im, and 77 is the cir- 
rus fraction of Xir. In Eq. (^, / and e depend on the 
dust-to-gas ratio]]. According to Hirashita (1999a, b), the 
dust-to-gas ratio is expressed as a function of metallicity 
(see also Dwek |1998t Lisenfeld & Ferrara p^98| ) . Therefore, 



we find easily that the relation between Ljr and SFR de- 
pends on metallicity via / and e. 



^ In this paper, we focus our main attention on only / and 
e, although we should keep in mind that the (1 — rj) may de- 
pend on the dust-to-gas ratio. The effect of (1 — r;) on SFR is 
discussed in Appendix A. 
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In this paper, we quantitatively examine the impor- 
tance of the metallicity to our conversion law, the conclu- 
sion being that the metallicity dependence is not always 
negligible when we compare SFRs of very young galax- 
ies and those of present galaxies. First, we examine the 



formula by IHKOO in Sect, 
to-gas ratio. Then, in Sect. 



focusing on the effect dust- 
J, we present the metallicity 
dependence of the IHKOO 's formula. In Sect. ^, we extend 
our discussion to the cosmic star formation history. Fi- 
nally, we summarize our conclusions in Sect. ^ . 

2. Dependence on dust-to-gas ratio 

In this section, we derive the metallicity dependence of the 
conversion formula (Eq. @)). It is convenient to define the 
factor CiR, as 



SFR 

CiR 



3.3 X IQ- 



10, 



0.4- 0.2/-f 0.6e 



yr 



(2) 
(3) 



We focus on the dependence of Cir on metallicity. The 
dependence of Cir on dust-to-gas ratio is incl uded through 
/ and e as will be described in Sects. 2.1 and P^. We 



will check our formula observationally in Sect. 2.6. Then, 
using the relation between the dust-to-gas ratio and the 
metallicity, we will obtain the metallicity dependence of 
Cir in the next section. We assume Case B (Baker & 
Menzel 1938), where the optical depths for the Lyman 
series are large enough. 

2.1. Dependence of f on dust-to-gas ratio 

The dependence of / ( Eq. ( ^) on dust-to-gas ratio is 
obtained from Spitzer ( 1978| , hereafter S78). We define 
rs,d and rs,H as the optical depths of the dust and of 
the neutral hydrogen atoms, respectively, for the Lyman- 
continuum photons over a path length equal to the 
Stromgren radius, rs (see S78 or Eq. (^ for the defini- 
tion of rs). Since the radius of a dust-free H ii region is 
estimated to be rs, it is useful to normalize the length 
scale by rs. We assume spherical symmetry and spatial 
uniformity of H ii regions in this paper for the simplicity. 

First, we calculate Ts,d/Ts,H for an H II region. Here, 
both Ts,d and rs,H are approximated with the optical 
thickness at the Lyman limit (912 A). Then we obtain 



TS,d 
TS,H 



TS.d 

TiHsrs 



12- 



B-V 



iV„s 



(4) 



where the dust extinction at 912 A, A912, is taken to 
be IHEb-v niag according to the Galactic extinction 
curve, and rs,d = ^912/(2.6 logj^Q e). Moreover, nn and 
s = 6.30 X 10"^^ cm^ denote the number density of the 
hydrogen and the absorption cross section for a hydrogen 
atom in the n = 1 level (Eq. (5.6) of S78), respectively, 
and the column density A^h is defined as iVn = riH^s. 



If the physical properties of grains such as the extinc- 
tion curve are unchanged, Eb-v/N-h is proportional to 
the dust-to-gas mass ratio V. According to S78, V = 
6x10-^ when Nu/Eb-v = 5.9x10^1 mag'^ cm'^ (Sects. 
7.2 and 7.3 of S78). Here, we apply these values to H 11 
regions. That is, we assume that the dust-to-gas ratio in 
the H II region is the same as that in the mean value in 
the interstellar space. Then, Eq. (H) reduces to 



1 



TS.d 

rs H 3100 V 6 x 10 



V 



(5) 



This is the same as Eq. (5.23) in S78 but an explicit ex- 
pression of the dependence on T). 

Next, we calculate rs,H = "■hs^'s- The Stromgren ra- 
dius is determined as 



rs = 1.4 



V 1048 



1/3 



-2/3 



102 cm- 



(6) 



where represents the number of ionizing photons emit- 
ted from central stars per second. In this equation, we have 
assumed that the temperature of the H 11 region is 8000 
K (i.e., the recombination coefficient to the n = 2 level is 
a(2) = 3.09 x 10^13 21^3 g-i. ^5^4) ^^y^^ 5 2 

of S78) and that — Jin, where n^, is the number den- 
sity of electrons, is satisfied within the Stromgren radius. 
Using the estimation of the Stromgren radius above, we 
obtain 



Ts,H = 2.7x 103( 



102 cm-3 



L/3 / iV, 



V1048 S- 



1/3 



(7) 



Combining this with Eq, 
V 



Ts.d is estimated as 



T'S.d 



0.87 



6 x 10- 



)( 



102 cm 3 



1/3 



V 1048 S- 



1/3 



(8) 



We note that rs, d becomes large as iVu increases. This is 
because the probability that the dust grains absorb the 
photons inside an H 11 region increases as the size of the 
region becomes larger. 

Next, we estimate the fraction of the ionizing photons 
absorbed by dust grains. Due to the grain absorption, the 
size of an H 11 region is smaller than rs. In other words, 
if we define as the ionization radius normalized by rs, 
y\ < \. (Without dust grains, y\ — 1.) A useful relation 
between rs^d and y\ is given in Table 5.4 of S78, where y\ 
is determined from the following expression (Eq. (5.29) of 
S78): 



V e 



dy = l 



(9) 



Using j/i, the fraction of the ionizing photons absorbed by 
hydrogen, /, is expressed as 



/ = 2/f 



(10) 
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This is the same as / in Eq. (||). In Fig. we show / as 
a function of rs^ d- 

When we consider the dependence of / on the dust-to- 
gas ratio, a large ambiguity exists: The number of ioniz- 
ing photons per H ii region is unknown, since the typical 
number and mass function of OB stars in an H li region is 
difficult to determine exactly. Fortunately, for the purpose 
of finding the dependence on 7), this is resolved by cali- 
brating the "Galactic" / with the value of Orion Nebula 
(Petrosian et al. 1972). According to them, / = 0.26 in the 
Nebula. They also commented that the value explains the 
IR emission from H ii regions. Adopting / — 0.2 6 as th e 
typical value of the Galaxy (see also Inoue et al. 2000b| ), 
we obtain Ts,d — 2.7 for the typical Galactic H ii regions 
(we consider 2? ~ 6 x 10^"^ for such regions) from Fig. |l|a. 
Hence we write 



rs,d -2.7 



V 



6 X 10-3 



(11) 



This is consistent with Eq. (g) if we assume riu — 
10^ cm~3 aj^j ^ 3 X 10''^ s-i. Thus, the net effects 
of mass function and number of OB stars are included in 
the value of Nu by adopting Eq. (p]). This simplicity is 
meaningful for our motivation to find the dependence of 
CiR on v. Combining Eq. ( ^t|) with Eq. (^, we obtain / 
as a function of V as shown in Fig. |l|b. 
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2.2. Dependence of e on dust-to-gas ratio 
In IHKOO, e in Eq. (||) is defined by 



e = l-e- 



(12) 



where Tnonion IS the mean optical depth of dust for non- 
ionizing photons. That is, e represents the efficiency of 
t he dus t absorption of nonionizing photons. Inoue et al. 
( 2000c ) estimated e = 0.9 (rnonion — 2.3) from the av- 



eraged visual extinction of Usui et al. ( |l998|) 's sample 
{Ay — 1 mag) and the Galactic extinction curve be- 
tween 1000 A and 4000 A by Savage & Mathis ( [19791 ). 
This wavelength range is fit for our purpose, since most 
of the nonionizing photons from OB stars are emitted in 
the wavelength range much shorter than 4000 A. The ef- 
ficient extinction in this short wavelength makes e large 
(i.e., nearly unity). 

It is obvious that Tnonion depends on the dust-to-gas 
ratio. If the column density of gas contributing to the ab- 
sorption of nonionizing photons is fixed, Tnonion is propor- 
tional to the dust-to-gas ratio. Since we are interested in 
the dependence of Cir on V, we simply adopt Tnonion oc V. 
Here, we determine the numerical value of Tnonion as 



= 2.3 



V 



6 X 10-3 



(13) 



so that Tn 



becomes 2.3 for the Galactic dust-to-gas 



ratio. By combining Eqs. ( |l^ ) and (p^), we obtain e as a 
function of T> as shown in Fig. 



Fig. 1. a Fraction of the ionizing photons absorbed by 
hydrogen, /, as a function of Ts.d (optical depth of dust 
over a distance equal to the Stromgren radius), b Fraction 
of the ionizing photons absorbed by hydrogen, /, as a 
function of the dust-to-gas ratio, T). 



2.3. Cirrus Fraction 

The fraction of the cirrus component, 77, remains to be de- 
termined. In this paper, an empirical value of r] is simply 
adopted. According to Lonsdale Persson & Helou (1987), 
rj ~ 0.5-0.7 for their sample spiral galaxies. We use 77 



0.5 

as adopted in IHKOO (the value averaged for Usui et al. 
1998's sample) for the current estimate for our SFR or 
Cir. As a first step, we adopt a constant ry, since a the- 
oretically reasonable determination of it is very difficult. 
The variation of the cirrus fraction by the change of T) 
will be briefly considered in Appendix A. We note that 
the range of Cir becomes larger if we consider the metal- 
licity dependence of the cirrus fraction (Appendix A). 

2.4. Absorption of Lyman-a photons by dust 

Since IHKOO 's derivation of Eq. (|^) is based on Case B, 
Lya photons are assumed to be easily trapped in an H 11 
region (Osterbrock 1989 ). Thus, during the resonant scat- 
terings in an H II region, all the Lya photons are assumed 
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Fig. 2. Fraction of the nonionizing photons absorbed by 
dust, e, as a function of the dust-to-gas ratio, T). 



to be absorbed by grains in IHKOO, which this paper is 
based on. 

If the dust-to-gas ratio is significantly smaller than the 
Galactic value, the dust grains might not efficiently absorb 
the Lya photons in H ii regions. However, H i envelopes 
on a galactic scale generally exists around H ii regions. A 
Lya photon is absorbed in the H i envelo pes so on after 
it escapes from an H ii region (Osterbrock 1961 ). Hence, 
even if there is only a small amount of dust, we expect 
a sufficient chance for dust to absorb the Lya photons in 
the H I envelope. We check this point below. 

First of all, we define a path length of the Lya photons 
as /path- It is estimated as being Zpath ^ ^l^Ji^ya, where 
TLya means the optical depth for Lya photons and /Lya 
is the mean free path of the Lya photons. The square to 
TLya means that the resonant scattering is assumed to be 
a random- walk process of photons. Next, we define an op- 
tical depth of dust grains for Lya photons over the length 
of I 

path Tdust - 

It is estimated as rdust ~ Tra^Zpath^'dust, 
where a is the size of grains and n^-ast is their number den- 
sity. Here, we will discuss whether the Lya photons have 
a chance to escape from galaxies. We are interested in a 
star- forming region surrounded by H i gas envelopes. The 
scale length of the H i envelope, L, may be estimated to 
be about 100 pc, which corresponds to the thickness of the 
disk of spiral galaxies. In such a case, TLya is estimated to 
be about L/^Lya- Then, we find rdust ^ Tra^ndust^^/^Lya- 
Using I?, it is expressed as 



T'dust 



Lya^dust 



(14) 



where mp is the proton mass and TTidust is the dust mass 
for an assumed spherical dust with mass density of 3 g 
cm^^. Adopting n-n = 0.1 cm~^ for a diffuse H i medium, 
a = 0.1 /zm, 2? = 6 X 10~^ (0.01 times smaller than the 
Galactic value), L = 100 pc, and the resonant-scattering 
cross section of Lya photons 2 x 10~^^ cm^ (i.e., ^Lya = 5 x 
lO^'^ cm), we find rdust ^ 5x 10^. Thus, rdust is much larger 



than unity. This means that most of the Lya photons are 
absorbed even if the dust-to-gas ratio is as small as 0.01 
times the Galactic value. 



2.5. Dependence of Cm, on dust-to-gas ratio 

In the above subsections, we have expressed / and e as 
a function of the dust-to-gas ratio (V) , while rj is treated 
as a constant (77 = 0.5). We also assume a typical star- 
forming region whose mean density of gas and production 
rate of ionizing photons are about 10^ cm^'^ and 3.0 x 10^^ 
s~^, respectively. Then, we can express Cir defined in Eq. 
(||) as a function of the dust-to-gas ratio. In Fig. ^, we 
present Cir as a function of V. From this figure, we find 
that the coefficient of the conversion from IR light to SFR 
becomes 4-5 times smaller for 2? ~ 6 x 10~^ (the Galactic 
value) than that for P ^ 6 x 10~^ (0.01 times the Galactic 
value). 

Cir is also changed if we adopt a different extinction 
curve. In the previous discussions, we have adopted the 
extinction curve of the Galaxy. The current paragraph 
examines how the relation between Cir and V changes 
when we adopt the extinction curve of the Small Mag- 
ellanic Cloud (SMC) with the other quantities adopted 
in this paper unchanged. In the SMC-type extinction law, 
the ratio of the optical depth in Eq. (||) becomes two times 
larger, because the UV extinction is enhanced in the SMC 
extinction. This indicates that the same estimation per- 
formed in Sect, ^is possible if we make the normalization 
of V half. Thus, we can find a rough dependence of / on 
V from Fig. |l]b by doubling V (i.e., the line in Fig. is 
moved to the left by 0.3 dex). Since V of SMC is about 
5 X 10-4 (Issa et al. |1990| ), we find / ~ 0.75 for the SMC 
extinction case, while / ~ 0.85 if we adopt the extinction 
law of the Galaxy. We can apply the same scaling to e. 
For I? ~ 5 X 10"'', e ~ 0.4 for the SMC extinction case, 
while e ~ 0.2 for the Galactic extinction case. As a result, 



at X> ~ 5 X 10-4, Cir ~ 3.4 x 10" 



-10 



M, 



yr 



for the 



SMC extinction, while Cir ~ 4.7 x 10"i° Mq yr"! Lq^ 
for the Galactic extinction. Thus, we should be aware of 
the effect of extinction law on our conversion formula, but 
the effect is small (a factor of 1.4). If we remember that 
the uncertainty owing to the IMF is a factor of 2 (e.g., 
IHKOO), we find that the dependence of Cir on the the 
extinction curve is not very important. We note that we 
can estimate Cir for the SMC-type extinction by moving 
the line in Fig. || to the left by 0.3 dex (a factor of 2). 

2.6. Observational check 

Here, we examine whether the values of / and e that we 
adopted for the Galactic value (i.e., / = 0.26 and e = 0.9) 
are consistent with the properties of actual H 11 regions, 
by using the ratio between the ionizing-photon luminosity 
and the IR luminosity. This r atio w as extensively studied 
in 1970s (e.g.. Harper & Low |l97l|) . 
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Fig. 3. CiR (conversion coefficient from the IR luminosity 
of dust to the star formation rate) as a function of the 
dust-to-gas ratio, V. 



We start with the relation by Petrosian et al. (1972) 



as 



L 



IR 



i(Lya) + (l-/)(HioniVu + ei 11 



(15) 



where {hv)ioa is the mean energy of an ionizing photon, 
and LiR, L(Lya), and inonion are luminosities of IR dust 
emission, Lya, and nonionizing photons, respectively. We 
note that the meanings of / and e in Eq. ( p^ are the same 
as those in Eq. (||). IHKOO derived Lnonion = l-5(/iJ^)ion-^u 
by adopting the Salpeter's initial mass function (Salpeter 
1955) with the stellar mass range of 0.1-100 Mq and the 



mass-luminosity relation of stars (Table 3.13 of Binney & 
Merrificld |1998| ). Then, according to Mezger ( |197S| ), we 
define the IR excess (IRE) as 



IRE = 



L 



IR 



= a+(l-/ + 1.5e) 



{hv)ion 
hvLyaf 



(16) 



where lyhya is the frequency of the Lya radiation, and a 
is the fraction of Lya photons that reach the 2p state and 
go down to Is (i.e., L(Lya) — ahi^i^yafNu^- According 
to S78, a ~ 0.67. Adopting / = 0.26 and e = 0.9 for the 
Galactic values (Sects. and 2^), we obtain IRE =11, 
where we have assumed that (ft.i^)ion is equal to hv at the 
Lyman limit (912 A). According to H arper & Low (1971), 
5 < IRE < 10 (see also Aannestad |1978| ; Mezger |l 9 78 ; 
Maihara et al. |l98lD for Galactic H ii regions. Fig. 7.3 
in S78 showed that the IR luminosity of H ii regions is 
larger than predicted from Lya luminosity by an order 
of magnitude (i.e., IRE ~ 10). Thus, / and e adopted in 
this paper are almost consistent with the properties of the 
Galactic H ii regions within the scatter of observed IRE. 

Another test of our model is possible for the case of H ii 
regions in the Large Magellanic Cloud (LMC) by using the 
observational results of DeGioia-Eastwood ( |1992D . From 

^ The other photons decay to the ground state with a two- 
photon process. 



Table 3 of DeGioia-Eastwood (1992), we obtain system- 
atically lower IRE for the LMC H ii regions; IRE ranges 
from 1 to 4, and the median is 1.6. This can be naturally 
explained by our model with a lower value of the dust-to- 
gas ratio in the LMC. If we adopt T) — 1.2 x 10~^ for the 
LMC (one-fifth of that of the galaxy) according to Issa et 
al. (1990), we obtain / ~ 0.7 and e ~ 0.4. Then we obtain 
IRE ~ 2.4 for the LMC. This value of IRE for the LMC 
lies in the range of the observed IRE. Thus, we consider 
Eqs. ( pi]) and (^3|) to be applicable and useful in spite of 
simplification in their derivation. Therefore, we discuss the 
metallicity dependence of Cm and its effect on the cosmic 
star formation history by using the relation between Cir, 
and V as shown in Fig. ^ in the following sections. 

3. Metallicity dependence 

In order to obtain the metallicity dependence of the con- 
version formula (Eq. (|l|)), which depends on the dust-to- 
gas ratio, 2?, as shown in the previous section, we must re- 
late V and the metallicity. Indee d, V and metallicity are 
posit ively correlated (Issa et al . 1990; Schmidt & Boiler 
|l993 ; Lisenfeld & Ferrara 1998). Here, we adopt the lat- 
est relation proposed by Hirashita (1999a, b). He has con- 
structed a new evolution model of dust in the galactic envi- 
ronment, emphasizing the importance that the grains can 
grow via the accretion of the metal elements in the cool 
and neutral components of interstellar medium (ISM; see 
also Dwek 199^ . Then, his relation between V and metal- 
licity has explained the observational relation for both the 
giant and dwarf galaxies. Here, we adopt the relation as 
shown in Fig. ^ (solid line), since this seems to repro- 
duce the relation for dwarf galaxies and spiral galaxies^ 
The metallicity is represented by [0/H]. We note that 
[0/H] = X means that the abundance of oxygen is 10^ 
times the solar value (the solar oxygen abundance is as- 
sumed to be 0.013 in mass; e.g., Whittet 1992). In Fig. 
^, observational data are also shown: The observed re- 
lations for nearby spiral galaxies are presented with the 



square (Issa et al. 199C), and the area marked with dwarfs 



represents a typical locus of dwarf irregular galaxi es an d 
blue compact dwarf galaxies (Lisenfeld & Ferrara 1998 ). 
We hereafter adopt the solid line in Fig. ^ as the relation 
between the dust-to-gas ratio and the metallicity. 

Once we accept the relation in Fig. we relate Cir and 
[O /H] by using the relation between the dust-to-gas ratio 
and Cir (Fig. ||). The relation between Cir and [0/H] is 
presented in Fig. ||. We see that if [O /H] in the ISM evolves 
from —2 to via chemical evolution, the coefficient of the 



As for the parameters in Hirashita (1999a, b), we choose 
the values as /in,o = 0.1 and /9acc = 2/3sn = 10, where /in,o 
is the dust mass fraction in the material injected from stars, 
and /3acc and /3sn are defined as gas consumption timescale 
(gas mass divided by star formation rate) normalized by dust 
growth timescale in clouds and by dust destruction timescale in 
supernova shocks, respectively (see Hirashita 1999a for details). 
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Table 1. Metallicity and Cir, as functions of z. 
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Fig. 4. Relation between the dust-to-gas ratio T) and oxy- 
gen abundance [0/H]. The solid line represents the well- 
fitting model by Hirashita (1999a, b). The observed re- 
lations for nearby spiral galaxies are presented with the 
square (Issa et al. 1990). The area marked with dwarfs 
represents a typical locus of dwarf irregular galaxies and 
blue compact dwarf galaxies (Lisenfeld & Ferrara 1998). 
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Fig. 5. Cir (conversion coefficient from the IR luminosity 
of dust to the star formation rate) as a function of the 
oxygen abundance, [0/H]. 



conversion from IR light to SFR becomes about 4-5 times 
smaller, again. Thus, when we would like to determine the 
SFR precisely within a factor of 4-5, we should not neglect 
the effect of metallicity. 

4. Comments on the cosmic star formation 
history 

The cosmological evolution of SFR is derived from the co- 
moving density of galactic light. For example, Madau et 
al. (1996) applied the conversion formula from UV light 
to the SFR and showed that the SFR as a function of the 
redshift z seems to have a peak at z ^ 1-2. The cosmic 
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has been revised and discussed 
1999D , commonly re- 



SFR as a function of z 
by many authors (e.g., Steidel et al 
ferred to as the Madau plot. Takeuchi et al. (2000) applied 
the formula by IHKGG to the determination of the cosmic 
star formation history from the number-count data of IR 
galaxies. 

The cosmic SFR is also determined from the dust emis- 
sion at the sub-millimeter (sub-mm) observational wave- 
length (Hughes et al. |199^ ; Barger et al. |1998D . We can 
convert the sub-mm luminosity to the SFR by applying 
the formula proposed in IHKOO. If we take into account the 
chemical evolution of galaxies on a cosmological timescale 
(e.g., Pei & Fall 1995| ), however, we must examine the de- 
pendence of the conversion formula on metallicity as seen 
in the previous sections. 

According to such a motivation, we apply the results 
obtained in the previous sections, especially Fig. || and its 
related discussions, to the cosmic star formation history. 
We consider the evolution of our conversion law along the 
metal enrichment with the aid of previous research that 
determined the metallicity as a function of z. We adopt Pei 



et al. (1999) as a recent modeling of the cosmic chemical 
evolution. 

Pei et al. (1999) modeled the cosmic star formation 
history and chemical evolution and calculated the evolu- 
tion of dust amount. Considering the absorption and re- 
processing of light by dust, they determined the cosmic 
chemical evolution after the calibration with the observed 
comoving density of light. The resultant metallicity evo- 
lution as a function of z is shown in the second column of 
Tab le p]. Here we adopt the solid line of Fig. 8 in Pei et 
al. ( 1999| ). Though their treatment of the evolution of the 
dust-to-gas ratio is not just the same as our treatment, for 
readers' qualitative understanding, we present a relation 
between the metallicity and the redshift z from their Fig. 
8. 

We also present Cir in the third column of Table |^ 
calculated from the relation between metallicity (Z) and 
Cir as shown in Fig. ^. Here we assume that XogiZjZQ) = 
[O /H] , where Zq is the solar metallicity and Z is measured 
in the galactic scale. From Table |^, we see that Cir is 4-5 
times larger at 2 3 than that at z ~ 0. Thus, we should 
carefully consider the chemical evolution of the galaxies 
if we determine the cosmic star formation history from 
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the dust emission within a factor of 4-5. Deriving the cos- 
mic star formation history from IR-sub-mm observation 
will be made possible by the future observational projects 
(e.g., Takeuchi et al. 1999| ). Our formulation, which takes 



into account the metallicity dependence, will be useful in 
determining the cosmic star formation history. 

We should mention that Pel et al. ( |1999|) treated aver- 
aged quantities for each redshift and focused only on the 
redshift dependence. Thus, our Table |^ must be applied 
to the data averaged for each z. In other words, we should 
not apply them to each individual galaxy at a certain z. 
For each galaxy, Fig. || should be used instead after their 
metallicity is known. Finally, we should keep in mind the 
possibility that far-IR (FIR) or sub-mm sample of galax- 
ies is biased to metal-enriched system because metal-poor 
(i.e., dust-poor) galaxies are not easily observed in these 
wavelengths. 



5. Conclusion 

Based on IHKOO's formulation, we consider the factor 
CiR in the conversion formula between the IR luminosity 
and the SFR (Eqs. (§ and (|)). The factor Cm becomes 
1.7 X 10-1°, 5.6 X 10-1°, and 7.9 x IQ-i" Mq yr"! L'^ 
for the metallicity of 1, 0.1, 0.01 times the solar value, 
respectively. Thus, Cm differs by a factor of 4-5 in the 
range. Importantly, applicability of our formula is obser- 
vationally confirmed by IREs of the Galaxy and the LMC. 
Applying our result to the cosmic star formation history, 
we have found that Cir may be about 4-5 times larger 
at z ^ 3 than that at z ~ 0. This factor of 4-5 is larger 
than the uncertainty caused by the initial mass function 
of stars (a factor of 2). Thus, we should carefully consider 
the chemical evolution of galaxies if we determine the cos- 
mic star formation history from the dust emission within 
a factor of 4-5. 
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Appendix A: dependence of cirrus fraction on 
metallicity 

In the main body of this paper, we have assumed that the 
cirrus fraction of the IR luminosity, 77, is 0.5. Here, we ex- 
amine how CiR, is affected owing to the change of 77. There 
is a large uncertainty about observational estimate of ry. 
But, theoretically, we consider some simple cases where 
the change in 77 cannot be neglected in our conversion for- 
mula. 




,0 -1, 



LO -0.5 
[0/H] 



0.0 0.5 



Fig. A.l. CiR (conversion coefficient from the IR dust lu- 
minosity to the star formation rate) as a function of the 
oxygen abundance, [O/H], for starburst galaxies. The cir- 
rus fraction 77 is assumed to be 0. 



The assumption of 77 = 0.5 breaks if we are interested 
in starburst galaxies. When we examine a sample of star- 
burst galaxies, indeed, it is reasonable to assume 77 ~ 0. 
This is because young stars dominate the radiation field 
that heats the dust, and the optical depth of dust is so 
large that almost all of the bolometric luminosity is emit- 

plsbj ). 



ted in the IR (Soifer et al. [1987t Kennicutt |1998bD . For 
the starburst galaxies, thus, the dependence of Cir on the 
metalHcity is obtained by putting 77 = into Eq. (||), while 
the dependence of / and e on the metallicity is the same 
as described in Sect. Isl. The result is shown in Fig. |A.l 



The result is shown in Fig 
We see that Cir in this figure is larger by a factor of 2 
than that in Fig. ||. 

Next, we consider a possible change in 77 of spiral galax- 
ies as a function of metallicity. The cirrus fraction of spiral 
galaxies is not negligible at the present epoch (77 ~ 0.5 fo r 
nearby spiral galaxies; Lonsdale Persson & Helou 1987). 
In the beginning of a spiral galaxy's evolution, when only 
a little metal is produced, 77 = 0.0 may be reasonable. 
Thus, the first approximation for the dependence of 77 on 
the dust-to-gas ratio is 



77 = 0.5 



V 



6 X 10- 



(A.l) 



If we adopt this relation, we obtain Fig. [A.2| as the rela- 
tion between Cir and metallicity. We note that our simple 



treatment as Eq. (A.l) is only applicable for [O/H] ^0.1, 
because 77 > 1 does not make sense. This means that we 
need a nonlinear modeling of the relation between rj and 
at the present epoch. We may need to consider a complex 
mode of evolution of ISM, whose nonlinearity causes inter- 



mittent star formation history (Kamaya & Takeuchi 1997; 
Hirashita & Kamaya |2000t Takeuchi & Hirashita pOOOD - 
Comparing Cir at [O/H] = -2 with that at [O/H] = 0, 
we see that Cir changes by an order of magnitude as the 
chemical enrichment proceeds. Thus, if we want to know 
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0.5 



[0/H] 

Fig. A. 2. Cm (conversion coefHcient from the IR dust lu- 
minosity to the star formation rate) as a function of the 
oxygen abundance, [0/H], for spiral galaxies. The cirrus 
fraction r] is assumed to be proportional to the dust-to-gas 
ratio. 



the realistic cosmic evolution of galaxies from Lm, the 
metallicity evolution of cirrus must also be examined. We 
will try this very difficult theme in the near future. 



Osterbrock D. E. 1989, Astrophysics of Gaseous Nebulae and 
Active Galactic Nuclei, University Science Books, Califor- 
nia 

Pei Y. C., Fall S. M., 1995, ApJ, 454, 69 
Pei Y. C., Fall S. M., Hauser M. G., 1999, ApJ, 522, 604 
Petrosian V., Silk J., Field G. B., 1972, ApJ, 177, L69 
Salpctcr E. E., 1955, ApJ, 121, 161 
Savage B. D., Mathis J. S. 1979, ARA&A, 17, 73 
Schmidt K.-H., Boiler T., 1993, Astron. Nachr., 314, 361 
Soifer B. T., Houck J. R., Neugebauer G., 1987, ARA&A, 25, 
187 

Spitzer L. Jr., 1978, Physical Processes in the Interstellar 

Medium, Wiley, New York (S78) 
Steidel C. C., Adelberger K. L., Giavalisco M., Dickinson M., 

Pettini M. 1999, ApJ, 519, 1 
Takagi T., Arimoto N., Vanscvicius V., 1999, ApJ, 523, 107 
Takeuchi T. T., Hirashita, H., 2000, ApJ, 540, 217 
Takeuchi T. T., Hirashita H., Ohta K., Hattori T. G., Ishii T. 

T., Shibai H., 1999, PASP, 111, 288 
Takeuchi T. T., Ishii T., Hirashita H., Yoshikawa K., Mat- 

suhara H., Kawara K, Okuda H., 2000, PASJ (in press) 
Tinsley B. M., 1972, A&A, 20, 383 
Usui T., Saito M., Tomita A., 1998, AJ, 116, 2166 
Whittct D. C. B., 1992, Dust in the Galactic Environment, 

lOP, New York 



References 

Aannestad P. A., 1978, ApJ, 220, 538 

Baker J. G., Menzel D. H., 1938, ApJ, 88, 52 

Barger A. J., Cowic L. L., Sanders D. B., Fulton E., Taniguchi 

Y., Sato Y., Kawara K., Okuda H., 1998, Nature, 394, 248 
Binney J., Merrifield M., 1998, Galactic Astronomy, Princeton 

University Press, New Jersey 
Buat v., Xu C., 1996, A&A, 306, 61 
DeGioia- Eastwood K., 1992, ApJ, 397, 542 
Dwek E., 1998, ApJ, 501, 643 
Harper D. A., Low F. J., 1971, ApJ, 165, L9 
Hirashita H., 1999a, ApJ, 510, L99 
Hirashita H., 1999b, ApJ, 522, 220 
Hirashita H., Karnaya H., 2000, A J, 120, 728 
Hughes D., et al., 1998, Nature, 394, 241 
Inoue A. K., Hirashita H., Kamaya, H., 2000a, PASJ, 52, 539 

(IHKOO) 

Inoue A. K., Hirashita H., Kamaya, H., 2000b, ApJ (submit- 
ted) 

Inoue A. K., Hirashita H., Kamaya, H., 2000c, AJ (in press) 
Issa M. R., MacLaren I., Wolfendale A. W., 1990, A&A, 236, 
237 

Kamaya H., Takeuchi, T. T., 1997, PASJ, 49, 271 
Kennicutt R. C. Jr., 1998a, ARA&A, 36, 189 
Kennicutt R. C. Jr., 1998b, ApJ, 498, 541 
Lisenfeld U., Ferrara A. 1998, ApJ, 496, 145 
Lonsdale Porsson C. J., Holou G., 1987, ApJ, 314, 513 
Madau P., Ferguson H. C., Dickinson M. E., Giavalisco M., 

Steidel C. C., Fruchter A., 1996, MNRAS, 283, 1388 
Maihara T., Oda N., Shibai H., Okuda H., 1981, A&A, 97, 139 
Meurer G. R., Heckman T. M., Calzetti D., 1999, ApJ, 521, 64 
Mezgcr P. G., 1978, A&A, 70, 565 
Osterbrock D. E. 1961, ApJ, 135, 195 



